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INTRODUCTION 


Metal matrix composites reinforced with ceramic fibers have high strength and 
stiffness to density ratios and are attractive for many applications where weight is a 
concern. They are also expected to have good high temperature strength. Advantage can be 
taken of the dissimilar properties of fiber and matrix to create composites with superior 
properties: for example ceramic fibers provide high temperature strength and a ductile metal 
matrix provides energy dissipation. However, ceramics have a low coefficient of thermal 
expansion (CTE) and metals have a high CTE. The mismatch in CTE induces thermal 
stresses in the composite when subjected to temperature change. Time variations of 
temperature can reduce the fatigue life of metal matrix composite significantly [1]. Brittle 
intermetallic matrix materials, that are attractive because of low density, could crack during 
the fabrication of the composite at cool down after consolidation [2]. Such cracks could 
impair the performance of the composite dramatically. 

Such problems have focused interest on the feasibility of using coatings on the fibers 
that would shield the matrix from the loading caused by the stiff fibers when the composite 
is subjected to a temperature change. Numerical parametric studies have been performed for 
continuous fiber reinforced composites by using concentric cylinder models [3] and by 
assuming the fibers are arranged in a square arrays [4], The studies suggest that the layer 
should have a very low modulus and a CTE that is the average of the CTE's of the matrix 
and fiber in order to reduce the residual stresses in the matrix. However, the conclusions are 
based on limited numerical results and do not provide an understanding of the fundamental 
problem. It appears that the formulation of a simple model which includes all the essential 
features of the problem and can yield insight into the physics of the problem has been 
overlooked. The results of such a model are not expected to be exact but provide guidelines 
on where to focus the subsequent numerical studies. 
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A simple model is formulated for a continuous fiber reinforced metal matrix 
composite with an interface layer when subjected to temperature change. The essential 
properties of an interface layer that could reduce the residual stresses in the matrix are 
identified. The possibility of improving the low cycle fatigue properties of the composite is 
also addressed. 


MODEL 

The composite consists of long aligned fibers with a thin interface layer in a metal 
matrix. The thermal loading of the composite is axi— symmetric with respect to the fibers 
and does not stimulate interaction between fibers in the transverse plane. The mechanical 
behavior of the composite may then be assumed to be given by a concentric cylinder model, 
Fig. 1. The composite is not subjected to transverse loading and the outer surface of the 
compound cylinder is traction free. The fibers are long and the stress and strain distributions 
are constant in the z— direction except at the end regions which are not studied here. 
Ceramic fibers are usually four to five times stiffer than metal matrix materials. Hence it is 
assumed, in order to simplify the calculations, that 

E f 

£ — >> 1 ( 1 ) 
m 


and 


p- » 1 (2) 

fc / 

where the subscript f denotes fiber, m matrix, and / interface layer. The fiber volume 
fraction , c^ , is assumed to be of the same order as the matrix volume fraction so that 
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( 3 ) 
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The radii defining the cylinder model (Fig.l) are related to the fiber volume fraction by the 
relation 


[a/b] 2 = c f 


( 4 ) 


The cylinder is assumed to be stress free at a given temperature T ^ and is then 
subjected to a homogeneous temperature change to a lower temperature T 2 • The 
temperature change is defined as 


AT = Tj 



( 5 ) 


For the given assumptions the stiff fiber controls the thermal expansion in the axial 
direction of the cylinder model. Hence, the stress distribution in the matrix and interface 
layer is governed by a plane strain problem with a fixed inner radius at the fiber interface. 
The problem is formulated in terms of the differential CTE’s of matrix and fiber 


Aa = a - a f 
m m 1 


(6a) 


and of interface layer and fiber 


ACL[ - otj - a f (6b) 

The stress distribution in the matrix is given by the plane strain solution for a thick 
walled cylinder subjected to an unknown internal pressure p , a traction free outer surface, 
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and a temperature change AT . Superimposing the solutions given in [5] for an internal 
pressure and a temperature change gives the stress distributions in the matrix 


a .-pD*t_l 

m [b/a] — 1 


CT _ _ [b/r] + 1 


(7a) 


(7b) 


a = p 
zm v 


2 v 

m 

[b/a] 2 - 1 


+ E Aa AT 
m m 


(7c) 


The hoop strain is 

e (l>m ' §~ ” — l lb/r)2+ 1 + (M 2 - ') 

t m [b/a] - 1 l 1 V m -I 

— [1 + v JAa AT (8) 

m m v y 

The interface layer is assumed to be sufficiently thin for the variation of stress in the 
radial direction to be neglected. Because of the rigidity of the fiber the deformation of the 
layer is suppressed in the hoop and axial direction by the fiber so that 




(9a,b) 


The layer is subjected to a pressure in the radial direction resulting from the contact with 
the matrix 
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(10a) 


°r/ = "P 

The temperature change AT introduces stresses in the hoop and longitudinal directions due 
to the constraint from the fibers (9). The stresses in the hoop and longitudinal direction are 
equal and are the sum of the stresses caused by the pressure and the thermal expansion 

V = °z / ■ - r p + E i Aa i AT ( 10b - c> 


The stresses in the layer give rise to the radial strain 


e r r~ 



(ID 


Compatibility in the radial displacement at the boundary between the inner surface 


of the matrix and the interface layer at r = a , requires 


v a = V 


( 12 ) 


The pressure in the radial direction at the layer is now derived by substituting (4) 
and (8) in (12) and using (11) to give 


t Aaj 1 + v / 


[ 1 + v m ] ' nr r-nr 

m / 


[1 - 


vi 


1 + c f , v m 
T ^{ 


tf 


m 


T^ E m A(X m AT < 13) 
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STRESSES IN MATRIX 

The highest stresses in the matrix occur at the inner surface of the cylinder, r = a. 
The principal stresses are given by (7a-c) by letting r = a and using (4). The v. Mises 
effective stress 


o*-r? 


+ a 


<P 




a a 
r <p 


- a r a z " 


a a 
<p z 


(14) 


can be calculated from the principal stresses given in (7) as 


°m = 


n f 
P 1 - C 


fJ 


•[3/0* + 1 -4v m (l - v m )l - p E m Aa m AT 2[1 - 2vJ + 


1/2 

[E Aa AT] Z 
m m J 


(15) 


The stresses in the matrix are readily found by first calculating the interface layer 
pressure for a given temperature change AT using (13) and using this value in (7a— c) and 
(15) to determine the stresses. The term E m Aa m AT enters in all the expressions for stress 
and all results can be normalized with respect to this term. The stresses at the inner radius 
of the matrix are shown in Fig. 2 as a function of normalized pressure for a fiber volume 
fraction c^ = 0.4 in a Ti^Al matrix for which v m = 0.25 . It can be seen that pressure at 
the interface causes increased tensile stresses in the hoop and axial direction which could 
result in matrix cracking, which is initiated from defects at the fiber matrix interface. The 
effective stress is a minimum when the interface pressure that is close to zero. Consequently 
the pressure should be low in order to have good low cycle fatigue properties of the 
composite. 
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For any given composite a plot similar to Fig. 2 can be produced. The value of the 

critical stress for matrix cracking , a cm , sets an upper limit for the principal stresses 

a and a and defines bounds on the pressure p if matrix failure is to be avoided. 
’ zm ’ rm 

The allowable effective stress , a ym , in the matrix defines other bounds on the pressure. 
These bounds define a window for the interface pressure (Fig 2.) if the failure criteria are 
not to be violated. If the calculated pressure does not fall within the allowable window it 
must be adjusted to avoid matrix cracking or failure due to high effective stress. 


Inspection of (13) reveals that the pressure may be reduced in two ways by the 
introduction of interface layers. A compliant layer with low modulus would increase the 
denominator in (13) while a compensating layer with a high coefficient of thermal 
expansion would decrease the numerator in (13). An appreciation for the potential of the 
two effects can be obtained by evaluating the requirements to reduce the pressure by a 
factor two. To achieve this condition using a compliant layer means that the two terms in 
the denominator of (13) are equal and gives the result, 

E / t V 1 ~ C f } (1A) 

E^'an - v ; )(l + v m )l(l - v m )(l + c f ) + v m d - c f >J 


For the present composite, using = 0.3 , this requires that 


E ; 


0.17 i 
a 


In general the layer has to be thin in order to keep the density down of the composite, and a 
realistic value of t/a is 0.1. This requires that the modulus of the layer must be 
unrealistically low for the pressure to be reduced by a factor two. Hence it is deduced that a 
compliant layer is not an efficient way to reduce thermal stresses in the matrix. 
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Furthermore, the second term in the denominator of (13) is usually sufficiently small to be 
neglected in the calculations and the pressure is then linearly dependent on the normalized 
thermal expansion of the layer, Acy/Aa m a . This is demonstrated in Fig, 3 where it can be 
seen that the influence of E m t/E^a is small. 


A reduction of the pressure by a factor two with a compensating layer requires that the 
second term in the numerator of (13) is equal to half of the first term, to give. 


Aa , 1 - v, 1 + v 

/ _ a / m 

Aa t 1 + v, 2 

m / 


(17) 


For the present composite this gives the condition 


Aa , 

SoT = 0 34 7 
m 


Hence when t/a = 0.1 the differential CTE for the layer is Aa, = 3.4 Aa . For the case of 
SiC fibers (a f = 5 10“^ 1/C) in a Ti^Al matrix (a m = 10 10 -6 ) requires that ~ 20 10 -6 
if the pressure is to be reduced by a factor two. Handbook values [6] of metals with a 
reasonably high melting point and high CTE are: silver 25 10 -6 , copper 17 lO - ^, and 
manganese 22 10~ It is evident therefore that an interface layer of readily available 
materials with high CTE has the potential to substantially reduce thermal stresses in the 
matrix. 


The results given in Fig. 2 and Fig. 3 have been combined in Fig. 4 where the 
highest principal stress and effective stress at the inner radius of the matrix cylinder are 
plotted against the normalized thermal expansion of the layer. No interface layer 
corresponds to Aa^t/Aa m a = 0. It can be deduced that a layer with a lower CTE than the 
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matrix will increase the stresses in the matrix and a layer with a higher CTE than the matrix 
reduces the stresses in the matrix. The effective and principal stress have minima for 
different layer characteristics and the required thermal expansion of the layer, A<y , has to 
be a compromise as discussed earlier. The highest effective stress in the matrix, (15), 
governs the low cycle fatigue performance of the composite when it is subjected to thermal 
cycling. Shake down to linear elastic stress strain response is expected when the cyclic 
effective stress, based on a linear elastic calculation, is less than twice the yield stress. The 
shake down condition is plotted in Fig. 5 as a function of the normalized thermal expansion 
of the layer. It can be seen that an interface layer can double the allowable temperature 
range for elastic shake down in the matrix. 

STRESSES IN INTERFACE LAYER 

The tradeoff for the improved stress state in the matrix is that the thermal mismatch 
has to be taken up by the interface layer. The deformation of the layer is constrained in the 
hoop and longitudinal direction so that tensile stresses develop in these directions. The 
highest principal stress in the interface layer which is given by substituting (13) in (10b), is 
shown in Fig. 6 as a function of the dimensionless parameters defining the problem. The 
stress is strongly dependent on the layer modulus, E^/E m , and CTE , Aa^Aa m , but 
weakly dependent on the layer thickness t/a . This indicates that the stress is dominated by 
the constrained thermal expansion [ given by the second term in (10b)] and is weakly 
dependent on the interface pressure [ given by the first term in (10b)]. The product Acx^t 
governs the reduction of stress in the matrix (Fig. 4) and must have a required value to 
reduce the stresses in the matrix to an acceptable level. It can be deduced from Fig. 4 and 
Fig. 6 that it is more favorable to have a thick layer and a moderate high layer CTE than a 
thin layer and a high layer CTE, to meet the requirements of stress reduction in the matrix 
and avoidance of high tensile stresses in the layer. 
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The effective stress in the layer is given by inserting (10a— c) in (14) to give 



E^Aa^AT + (1 — 2vj) p 


( 18 ) 


Shake down to elastic response in the layer is expected when the cyclic part of the 
effective stress is less than twice the yield stress of the layer. The shake down limit has 
been plotted in a normalized form in Fig. 7 for different layer modulus and thickness. The 
stress in the layer is strongly dependent on the layer modulus and to ensure shake down it is 
desirable to have a layer with a high yield stress and a low modulus. 


THE EFFECT OF YIELDING OF THE INTERFACE LAYER 


The previous calculations indicate that the layer can be subjected to high stresses 
and may yield. An estimation of how yielding in the layer affects reduction of the stresses 
in the matrix can be obtained by examining the behavior of an elastic— perfectly plastic 
layer. In the analyses the layer is assumed to be in an initial elastic stress state which 
satisfies the yield condition (14). The cylinder is subsequently subjected to a further 
temperature change dT so that the radial pressure changes by dp . The loading is such 
that the layer remains yielded during the subsequent deformation requiring the effective 
stress to be constant. Differentiation of (14) and using the condition at yield, (10b,c), gives 
an incremental stress change in the layer which is hydrostatic. Thus, the stress increments 
are 


do z/ = d V = do rt ■ 


(19a,b,c) 
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Since the total strains of the layer in the <p and a directions are constrained to be zero, 
(9a,b), it follows that the plastic strains in these directions are the sum of the elastic and 

thermal part 


<( = de zr- io / dT 


1 -2v, 


dp 


(20a, b) 


The plastic deformation is incompressible requiring 

<21> 

and the total strain increment in the radial direction is the sum of the elastic, thermal, and 
plastic component and becomes 

de Pj = 1 ~ 2 ^ dp + Aa^dT + deP ; = 3Aa^dT - 3 E y — dp ( 22 ) 

The radial strain in the layer governs stress reduction in the matrix and an appreciation for 
the effect of the yielded layer can be determined by comparing the radial strain for an 
elastic layer, Eq. (11), with the radial strain for an elasto-plastic layer, Eq. (22). The 
contribution from the pressure p in both these equations is small and can be neglected. 
The strains are then only linearly dependent on the temperature change. By letting AT = 
-dT in Eq. (11) the ratio of the two strains is found to be 


de 

de 


PI 
r / 

r / 


3(1 - vp 
1 + v," 


( 22 ) 
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To estimate the effect of the elasto— plastic deformation, the second term in the numerator 
of (13) must be multiplied by the right hand side of (22) when calculating the pressure. For 
the present data , with = 0.3 , this implies that the thermal expansion of the layer 
becomes 1.5 times more effective in reducing the pressure. Hence the results of the elastic 
solutions given in Fig. 3—5 can be readily modified to include plastic deformation in the 
layer. 


CONCLUSIONS 

A compliant layer must to be unrealistically flexible to reduce thermal stresses in the 

matrix. 


A compensating layer with a high coefficient of thermal expansion (not 
unrealistically high) has the potential of reducing the thermal stresses in the matrix 
significantly. 

The hoop stress in the matrix close to the fiber can be reduced substantially but the 
axial stress in the matrix is less affected by a layer. This implies that compensating layers 
can be expected to be successful in preventing cracking in composites where predominantly 
radial cracking is observed in the matrix. 

The shake down range to linear elastic response in the matrix can be extended 
significantly with the addition of a compensating layer. 

The thermal mismatch between fiber and matrix is taken up by the layer which can 
cause high stresses in the layer. The required stress reduction in the matrix sets a 
requirement on the value of the thermal expansion of the layer, Acy . However, the stress 
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in the layer is proportional to Aot^ so that a high layer CTE can cause high stress in the 
layer. In order to have low stress in the layer it is more favorable to have a moderately high 
layer CTE and thick layer than a thin layer with very high layer CTE. A low elastic 
modulus of the layer reduces the stresses in the layer. The effect of plastic deformation of 
the layer on the stress reduction in the matrix is equivalent to increasing the CTE of the 
layer by 1.5. 

Thermal stresses build up in the layer because the required thermal expansion in the 
radial direction is different from the thermal expansions in the hoop and longitudinal 
directions that are governed by the fiber. The stresses in the layer could be minimized if a 
layer could be developed with an anisotropic thermal expansion such that the thermal 
expansion is close to the fiber CTE in the longitudinal and hoop directions and high in the 
radial direction. 

Residual stresses can be reduced by Compensating layers for other reinforcements 
shapes than fibers. Thermal expansion of reinforcements with lower aspect ratios is less 
anisotropic and can be better compensated for by a layer with an isotropic CTE. The thermal 
mismatch for a spherical inclusion can be compensated for fully and no residual stresses would 
develop in the matrix. 
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Fig. 2 Stresses at the inner surface of the matrix cylinder as a function of 

normalized pressure at the matrix Fiber interface. 


16 


p/E m Aa m AT 



A«|t/ Ao^nnQ 


Pig. 3 Relation between normalized pressure at matrix fiber interface and relative 

coefficients of thermal expansion and layer thickness. 
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Fig. 4 Maximum principal stress and effective stress at the inner surface of the 

matrix cylinder as a function of the relative thermal expansion of the layer. 
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